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Abstract 

The evolution of bacterial consortia was studied in six semi-solid rye sourdoughs during long-term bacl<slopping at different 
temperatures. Each rye sourdough was started spontaneously in a laboratory (dough yield 200), propagated at either 20°C 
or 30°C, and renewed daily at an inoculation rate of 1:10 for 56 days. The changes in bacterial diversity over time were 
followed by both DGGE coupled with partial 16S rRNA gene sequencing and pyrosequencing of bar-coded 16S rRNA gene 
amplicons. Four species from the genus Lactobacillus (brevis, crustorum, plantarum, and parallmentarius) were detected in 
different combinations in all sourdoughs after 56 propagation cycles. Facultative heterofermentative lactic acid bacteria 
dominated in sourdoughs fermented at 30"C, while both obligate and facultative heterofermentative LAB were found to 
dominate in sourdoughs fermented at 20°C. After 56 propagation cycles, Kazachstania unispora (formerly Saccharomyces 
unisporus) was identified as the only yeast species that dominated in sourdoughs fermented at 20°C, while different 
combinations of strains from four yeast species {Kazachstania unispora, Saccharomyces cerevisiae, Candida krusei and 
Candida glabrata) were detected in sourdoughs propagated at SO^C. The evolution of bacterial communities in sourdoughs 
fermented at the same temperature did not follow the same time course and changes in the composition of dominant and 
subdominant bacterial communities occurred even after six weeks of backslopping. 
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Introduction 

Sourdough is a mixture of flour and water fermented with a 
microbial community mainly consisting of lactic acid bacteria 
(LAB) and yeasts. LAB dominate the microbial community and 
are responsible for acid production while yeasts work to leaven the 
dough [1,2]. Lactic acid fermentation plays an important role in 
the production of rye bread by both decreasing the activity of ot- 
amylase and improving dough texture [3]. 

Lactobacillus brevis, Lactobacillus plantarum, and Lactobacillus 
sanfranciscensis are the LAB species most frequendy isolated from 
rye sourdoughs [2,4—5]. However, Lactobacillus amylovorus, 
Lactobacillus fermentum, Lactobacillus helveticus, Lactobacillus 
panis, Lactobacillus pontis, and Lactobacillus reuteri have also 
been identified as dominant species in rye sourdoughs [4,6-9]. 
Since rye flours have a generally higher extraction rate than wheat 
flours [5] rye sourdoughs are characterized by higher buffering 
capacity compared with wheat and spelt sourdoughs but also 
higher concentration of mannitol and amino acids [8,10]. Despite 
this, the same LAB species are able to dominate the microbial 



communities within mature sourdoughs made with different types 
of flour under otherwise identical conditions [8,10-1 1]. Recently, 
it has been shown that initial differences in the microbial 
communities found in spontaneously started rye and wheat 
sourdoughs decreases during backslopping propagations and that 
a common core microbiota is established [1 1]. 

The establishment and stability of microbial consortia in 
sourdoughs depends on the microbial communities within the 
raw materials, the chemical composition of the raw materials, and 
interactions between the microorganisms, together with fermen- 
tation parameters such as temperature, inoculum size, dough 
yield, and fermentation length [12-16]. Previous studies on the 
dynamics of microbial communities in spontaneously started 
sourdoughs have stated that the establishment of a stable 
consortium occurs through a three-stage evolution process within 
a few days during which one may observe the prevalence of 
sourdough-atypical LAB, sourdough-typical LAB, and highly 
adapted sourdough-typical LAB [2,8,10-1 1]. Most of these studies 
have been carried out with liquid sourdoughs fermented using 
continuous stirring and propagated over a short period of time 
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(maximum two weeks). Traditional Estonian rye sourdoughs can 
be classified as semi-solid because they have a relatively dense 
consistency. They are fermented at ambient temperature, which 

may vary between 18-27°C depending on the season. Because the 
water content (dough yield) and fermentation temperature are the 
main factors that affect the composition of the bacterial 
community [2], we evaluated the establishment of microbial 
consortia in spontaneously fermented rye sourdoughs with a low 
dough yield at two different temperatures (20 and 30°C) back- 
slopped daily for 56 days. The aims of our study are (i) to gain 
insight into the stability of bacteria communities after tludr initial 
establishment in mature sourdough, (ii) to determine the effect of 
fermentation temperature on the development of bacterial 
communities in sourdoughs, (iii) to assess the reproducibility of 
the development of microbial communities in sourdough when 
applying fermentation conditions with low dough yield. 

Materials and Methods 

Sourdough Fermentation and Sampling 

Sourdough fermentation was initiated in six 400 ml sterile 
Stomacher circulator bags (Seward Limited, England) by mixing 
150 g of rye flour and 150 g of sterile 0.5% NaCl solution. Rye 
flour (Type 1370, Tartu Mill, Estonia) from a single 50 kg bag was 
used during the entire experiment. Each sourdough was mixed for 
15 minutes at 100 rpm using a Stomacher 400 circulator (Seward 
Limited, England), the bags were sealed with tape, and positioned 
vertically into an incubator. Three sourdough batches were 
fermented at 20°C and another three at 30°C. After 24 hours of 
fermentation, each sourdough was thoroughly mixed in the 
circulator for 5 min at 100 rpm and renewed at an inoculation 
rate of 1:10 by mixing 30 g of sourdough, 135 g of sterile 0.5% 
NaCl solution, and 1 35 g of rye flour. In total, each sourdough 
experienced 56 backslopping cycles. The three sourdoughs 
fermented at 20°C are referred to as 20-1, 20-11, and 20-III while 
those fermented at 30°C are referred to as 30-1, 30-11, and SO-IIL 
Day 0 indicates the start of the experiment. 

Both the acidity, in units of pH, and total titratable acidity 
(TTA) were measured using a DL22 Food and Beverage Analyzer 
(Mettler-Toledo LLC, USA) at the end of each fermentation cycle 
prior to renewal. Following this, five grams of sourdough was 
homogenized with 50 ml of distilled water using a Polytron 
PT2100 homogenizer (Kinematica AG, Switzerland). 

Both bacterial plate counts and DGGE analysis were carried out 
on days 0, 1, 3, 5, 7, 10, 15, 21, 28, 35, 42, and 56 of tiie 
experiment. Pyrosequencing was performed on samples from days 
1, 3, 5, 7, 21, and 56. The sample from day 0 was taken from the 
rye flour and water mixture immediately after mixing. 

Isolation and Characterization of Lactic Acid Bacteria and 

Yeasts 

Five grams of sourdough were supplemented with sterile 0.85% 
NaCl solution up to a volume of 50 ml. The mixture was then 
homogenized by vortexing. Decimal dilutions were plated on to 
both sourdough bacteria (SDB) agar (maltose, 2.0%; yeast extract, 
1%; Tween 80, 0.03%, trypticase 0.6%; pH 5.6) and de Man, 
Rogosa and Sharpe (MRS) agar (Lab M Ltd, UK) with 100 |a.g/ml 
cycloheximide (Sigma-Aldrich, USA). 

The plates were incubated at the same temperature the 
sourdough was fermented at (20 or 30°C). Incubation was carried 
out for 48 h under anaerobic conditions (AnaeroGen, Oxoid). 
Colony forming units (CFU) were counted from the agar media 
using suitable dilutions. 



For each of the six sourdough samples collected on day 56, 20 
colonies were picked from the MRS and SDB agar plates (ten from 
each medium) for further analysis by rep-PCR. Colony picking 
was performed in succession from one sector of the plate. On day 
56 samples were also plated on Yeast Extract Peptone Dextrose 
(YPD) agar (dextrose, 2.0%; peptone 2%; yeast extract, 1%) with 
100 |J,g/ml chloramphenicol (Sigma-Aldrich, USA) and incubated 
at the same temperature the sample was fermented at (20 or 30°C). 
Ten colonies per sample were picked in succession from YPD agar 
plates and analyzed using RAPD-PCR. 

Extraction of Whole Genomic DNA 

Total DNA extraction was performed using 5 g of sourdough, 
which had been homogenized by vortexing with 45 ml of sterile 
physiological solution [10]. This suspension was then centrifuged 
at 4''C for 5 minutes at 1000 xg. The supernatant was collected 
and centrifuged at 4°C for 15 minutes at 5000 xg. Each extraction 
of whole DNA was performed using a GenElute Bacterial 
Genomic DNA Kit (Sigma-Aldrich., USA), according to the 
manufacturer's protocol. Total DNA was also extracted from the 
rye flour and water mixture immediately after mixing on day 0. 

Denaturing Gradient Gel Electrophoresis Analysis 

The V3 region of the 16S rRNA gene was amplified from the 
whole genomic DNA using universal primers F357-GC and 51 BR 
as described in [10]. DGGE was carried out using the 
INGENYphorU system (Ingeny International Bv., Netherlands) 
as described in Viiard et al. [9]. All clearly visible bands were cut 
from the gel and incubated in TE buffer (10 mM Tris pH 7.5 and 
1 mM EDTA pH 8.0) at 37°C for 1 h. Eluted DNA was 
reamplified using primers F357 and 518R and sequenced in a 
commercial facility (Estonian Biocentre, Tartu, Estonia). Nucleo- 
tide secjuences were analyzed using the BLASTn algorithm 
together with the GenBank database (National Center for 
Biotechnology Information, USA). 

Pyrosequencing of Bar-coded 16S rRNA Gene Amplicons 

Universal primers 8F and 357R were used for PGR amplifica- 
tion of the V1-V2 hypervariable regions of 16S rRNA genes [17]. 
The amplicon mixtures were pyr()sc'qu(;nced using a 454 GS 
FLX+ System (Roche 454 Life Sciences, USA) in a university 
facility (Centre for Biology of Integrated Systems, Estonia). 

The resulting pyrosequencing data was analyzed using the 
software package MOTHUR, version 1.27.0 [18]. Reads shorter 
than 150 bps were removed from the dataset and the PyroNoise 
algorithm was used to discard both homopolymer-derived and 
PGR errors. The remaining sequences were aligned to the SILVA 
reference 168 ribosomal RNA database [19]. Chimeric sequences 
were filtered using the UChime method by applying the 
'chimera.uchime' procedure in MOTHUR in de novo mode, which 
checks chimeras in each group separately. Operational Taxonom- 
ic Units (OTUs) were defined using the average neighbour 
clustering algorithm within MOTHUR with a 97% similarity 
threshold. Rarefaction curves and normalized OTU counts at 500 
sequences were calculated using the R software package "vegan" 
version 2.0-7. In addition, we calculated the rate of forming new 
OTL^s \\ licn one sequence is added to the set of 500 sequences. 
The closest match on the species level was found for each OTU 
using the BLASTn algorithm together with tlu- Gc-nBank database 
(National Center for Biotechnology Information, USA) with the 
parameters of 97% similarity and 90% coverage. The relative 
abundance of OTUs was calculated as the number of serjuences 
for each OTU divided by the total number of bacterial sequences 
obtained for each sourdough sample. 
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Statistical Analysis 

Plate count data was subjected to Z-tests to both compare the 
results obtained at 20 and 30°C and to compare results obtained 
using different media. We tested the hypotheses that the difference 
between the results obtained at different conditions come from a 
distribution with mean zero. The nuisance parameter of the z-test 
is the sum of the standard deviation of the differences between the 
values obtained at two conditions and the standard deviations of 
the values obtained at each condition. 

Acidity and TTA measurements obtained at 20 and 30°C were 
compared using a simple two sample Student's t-test for samples 
with equal variance. The values tested were considered to come 
from distributions with different mean values when the p-value 
was below 0.05. 

DNA Isolation from Colonies, Rep-PCR Fingerprinting and 

Partial Sequencing of rRna Genes 

Bacterial DNA was extracted from isolated colonies using 
Whatman indicating FTA MiniCards (GE Healthcare Ltd., UK) 
using a method provided by the manufacturer. Rep-PCR with 
primer (GTG)5 followed by agarose gel electrophoresis was 
performed as described by Viiard et al. [9] . Extraction of yeast 
DNA was carried out using a PureLink Genomic DNA mini Kit 
(Invitrogen, USA) using a method provided by the manufacturer. 
RAPD-PCR with an Ml 3 primer was performed according to 
Andrighetto et al. [20] . Each fingerprint type was calculated as the 
ratio of similar fingerprints to the number of colonies analyzed. 
One or two representatives from each PGR fingerprint group were 
subjected to Sanger serjuencing. The resulting 16S rRNA gene 
fragments were amplified using universal primers 27f-YM [21] 
and 16R1522 [22]. Yeasts were identified using standard protocols 
by amplifying the D1/D2 variable domains of the 26S rRNA gene 
with primer pair NLl and NL4 [23]. Amplified fragments were 
purified with a GeneJET PGR Purification kit (Fermentas, Vilnius, 
Lithuania). Sequencing of the fragments was conducted at a 
commercial facility (Estonian Biocenter, Estonia). The resulting 
gene sequences were compared with the GenBank database using 
the BLASTn algorithm (National Center for Biotechnology 
Information, USA). 

Evaluation of Carbohydrate Fermentation Profiles 

Carbohydrate fermentation profiles of selected LAB and yeast 
strains were determined with the identification kits API 50 CH and 
API 20 C AUX using methods provided by the manufacturer 

(bioMerieux, France). 

Results 

pH, Total Titratable Acidity, and Bacterial Plate Counts 

Spontaneous rye sourdoughs were started at 20°C and 30°C in 
three parallels. After 24 hours of fermentation at 20°C the average 
viable count of bacteria increased from 2.0x10' GFU/g in the 
raw flour to 5.7x10 while the first fermentation cycle at 30°G 
resulted in 7.2 x 10*^ CFU/g (Table 1). After the third backslopping 
cycle (day 3) the viable count of bacteria in all sourdoughs 
exceeded 10^ CFU/g. No significant (p<0.05) difference was 
found in the plate counts obtained from SDB or MRS media. 

Fermentation both increased the concentration of acids and 
lowered pH of the rye sourdoughs. For the fermentation series 
conducted at 30°C, the maximum TTA value reached (22.5±0.7) 
occurred on day 10 while the 20°G fermentation series displayed 
its maximum TTA value (19.5± 1.0) on days 1 1 and 12 (Figure 1). 
The maximum value for the 30°C series was larger, however, the 



difference in TTA between sourdoughs fermented at different 
temperatures diminished during later propagation cycles (Figure 1). 

During the first ten cycles, the pH in the sourdoughs fermented 
at 30°G decreased significandy more c:ompared with those 
fermented at 20°C (Figure 1). However, during the later stages 
of backslopping, the acidity did not significandy differ between 
sourdoughs. 

Dynamics of the Bacterial Communities Determined by 
DGGE of 16S rRNA Gene PCR Amplicons 

The highest diversity of species detected using DGGE analysis 
occurred after the first 24 hours of fermentation (Figure 2A). Over 
15 bands were visible, however, only few of the sequences obtained 
by cutting these hands were identified at a similarity level S97% 
(Table 2). Of those identified, three belong to Pantoea a^lomerans, 
one to a Gamma Proteobacteria species and one to a WeisselM species. 
All of the bands identified were observed at both fermentation 
temperatures. After the third backslopping, the DGGE patterns of 
all six sourdoughs were similar to each other and consisted of 
fragments specific to the LAB genera Weissella, Lactobacillus, and 
Pediococcus. On day 5 the DGGE patterns of sourdoughs fermented 
at 20°C had not drastically changed. However, from sourdoughs 
fermented at 30°C, the bands specific to Weissella sp. and 
Lactobacillus sp. had disappeared and another band, identified as 
Lactobacillus phmtamm, emerged (Figure 2A). 

On day 7, a band specific to L. plantarum was also detected in all 
sourdoughs fermented at 20°C together with another band 
identified as Pediococcus pentosaceus. In addition, a third new band 
identified as Lactobacillus cuwatus / graminis / sakei was found in batchs 
20-1 and 20-11 (Figure 2B), however, the DNA sequence obtained 
from this band did not provide sufficient information to 
discriminate between these three Lactobacillus species even at 
99% identity (Table 2). Differences were also found among 
sourdough batches fermented at 30°C. Fragments specific to 
Lactobacillus brevis and P. pentosaceus were seen only in sourdoughs 
30-1 and 30-III in addition to L. plantarum, which was found in all 
three batches. The DGGE pattern of 30-11 had two additional 
bands, one of which was identified as L. sakei (Figure 2B and 
Table 2). 

The diversity of the six bacterial communities continued to 
decrease during the second and third week of propagation 
(Figure 2B). On day 21, two major bands specific to the species 
L. plantamm and L. brevis were observed in sourdoughs 20-1 and 20- 
III while L. plantarum and Lactobacillus crustorum were detected in 
sourdough 20-11. In sourdoughs fermented at 30°C, only one 
single strong band corresponding to either Lactobacillus paralimentar- 
ius (30-1) or L. plantarum (30-11 and 30-111) was observed after cycle 
21 (Figure 2B). 

The bacterial composition of the sourdoughs was both stable 
and comparable between batches from day 21 to 35. However, 
further succession of species occurred after the fifth week of 
propagation (Figure 2C). On the final day of sampling, sourdoughs 
fermented at 20°C were comprised of either L. brevis and L. 
paralimentarius (sourdough 20-1), L. brevis, L. plantarum and L. 
crustorum (sourdough 20-11), or L. brevis and L. plantarum (sourdough 
20-III). In sourdoughs fermented at 30°C L. paralimentarius was 
detected in aU three sourdoughs either in combination with L. 
plantamm and L. crustorum (sourdough 30-1), together with L. 
plantamm (sourdough 30-11) or as a single dominant species 
(sourdough 30-III). 
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Cycle number 

Figure 1. pH and total titratable acidity of six rye sourdoughs propagated at 20°C (20-1, 20-11 and 20-111, sKiown in blue) and 30°C 
(30-1, 30-11 and 30-111, shown in red) during 56 days. 

doi:10.1371/journal.pone.0095449.g001 



Dynamics of the Bacterial Communities Determined by 
Pyrosequencing of Bar-coded 16S rRNA Gene Amplicons 

Pyrosequencing of bar-coded 16S rRNA gene amplicons was 
applied to overcome limitation of DGGE analysis and study in- 
depth the establishment of microbial consortia in spontaneously 
started rye sourdoughs. Widely used V1-V2 region specific 
primers were chosen to distinguish between different LAB species 
[24—26], After matching the barcodes and performing initial 
quality processing using PyroNoise, 48912 raw reads were found. 
After removing 1258 chimers and plant chloroplast related 
sequences, a total of 41819 high-quality partial 16S rRNA gene 
sequences longer than 150 bp were used in the analysis (Table SI). 
The majority of high-quality reads were in the length range 280- 
310 bp. The number of detected OTUs, normalized number of 
expected OTUs at 500 sequences, and the rate of new OTUs 
when one sequence is added to the sample at 500 sequences are 
provided in Table SI. The latter quantity indicates that the 
sequence coverage was sufficient for most of the samples, with the 
exception of four (0-20-1, 1-20-II, I-30-I, 5-30-III), which all had 
a lower number of reads (£500). During the first three 
propagation cycles, species diversity was higher in sourdoughs 
fermented at 30°C as determined by the normalized number of 
expected OTUs (Table SI). This number of expected OTUs per 
500 reads decKned during the continuous propagation of the 
sourdoughs to between 4 and 9 after 56 renewals (Table SI). 

Several bacterial species were found in the rye flour (Table S2), 
although the ratio of microbial DNA to rye DNA in this sample 
was low. As a consequence, only 60 high-quality bacterial reads 
were obtained. The main OTUs identified were from the genera 
Pantoea (33%) and Stmotrophomonas (15%) (Table S2). LAB species 
were represented by Lactobacillus iners (3.3%), Leuconostoc citreum 
(1.7%), and Wdssella cibaria (1.7%). 

After 24 hours of spontaneous flour fermentation at 20°C 
between 70-90% of all identified 16S rRNA gene amplicons were 
from the genera Enterobacter, and Pantoea (Table S2). Bacteria from 
the genera Wdssella, Leuconostoc, Lactococcus, and Lactobacillus 
together comprised up to 23% of the bacterial community. In 
sourdoughs fermented at 30°C, the majority of the bacterial 



community was formed by representatives of the genera Entero- 
bacter, Weissella, Lactococcus, and Leuconostoc. In contrast with the 
DGGE analysis, differences between the bacterial composition of 
sourdoughs fermented under the same conditions were already 
observed after the first fermentation cycle (Figure 3). 

In sourdoughs fermented at 20°C, enterobacteria were totally 
replaced by the LAB species Weissella cibaria, Weissella para- 
mesenteroides/confusa, Lactococcus lactis, Leuconostoc citreum, 
and Lactobacillus graminis/ sakei/ curvatus after the third renewal 
cycle on day 3. W. paramesenteroides/confusa and Lc. lactis 
formed the dominant population in sourdoughs 20-11 and 20-III. 
In batch 20-1, W. paramesenteroides/confusa and L. graminis/ 
sakei/ curvatus dominated. The proportion of L. graminis/ sakei/ 
curvatus increased in all sourdoughs fermented at 20°C after the 
fifth fermentation cycle. This species kept its dominant position in 
batches 20-1 and 20-11 up to day 8 and formed over 80% of the 
bacterial community. In contrast, L. plantarum/pentosus formed 
over 70% of the bacterial community in sourdough 20-III. After 
day 21, L. plantarum and L. brevis were the dominant species in 
sourdoughs 20-1 and 20-III, while in sourdough 20-11 L. 
crustorum and L. plantarum formed 85% and 13% of identified 
amplicons, respectively. Even after 56 days of fermentation the 
dominant community in all three sourdoughs fermented at 20°C 
diflered. In sourdough 20-1, Lactobacillus paralimentarius/kimchii 
(the latter is the synonym for L. paralimentarius [27]) and L. brevis 
dominated in the proportion 81:19. In sourdough 20-11 the 
dominant community was represented by L. crustorum and L. 
brevis in the proportion of 46:50. In sourdough 20-III L. 
plantarum and L. brevis were dominated in the proportion 34:65. 

In contrast with sourdoughs fermented at 20°C, enterobacteria 
were still present in low numbers within sourdoughs fermented at 
30°G after the third renewal cycle. Various combinations of W. 
paramesenteroides/confusa, Pediococcus pentosaceus, Leuc. citreum, and L, 
graminis/ sakei/ curvatus were found in sourdoughs 30-1, 30-11, and 
30-III. After the fifth renewal cycle, the dominant bacteria in 30-1 
and 30-11 were P. pentosaceus and L. plantarum/pentosus, respectively. 
In sourdough 30-III, L. plantarum, L. graminis /sakd/ curvatus, and 
Lactococcus lactis species were found in equal proportions (Figure 3). 
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Figure 2. DGGE analysis of six rye sourdoughs propagated at 20'C (20-1, 20-11 and 20-111) and 30°C (30-1, 30-11 and 30-111) during 56 
days. The arrows indicate specific bands for the following groups: (A) Weissella sp., (B) Lactobacillus plantarum, (C) Enterobacteriaceae, (D) 
Lactobacillus sp., (E) Lactobacillus sp., (F) Pediococcus sp., (G) Cereal chloroplast, (H) Pantoea agglomerans, (I) Lactobacillus brevis, (J) Lactobacillus 
paralimentarius, (K) Lactobacillus crustorum, (M) Lactobacillus curvatus/graminis/sakei, (N) Pediococcus pentosaceus, (O) Gamma Proteobacteria species. 
The bands marked by numbers were cut and sequenced. The closest matches are shown in the Table 2. A. DNA samples from days 0, 1, 3 and 5 of 
backslopping. B. DNA samples from days 7, 10, 15 and 21 of backslopping. C. DNA samples from days 28, 35, 42 and 56 of backslopping. 
doi:1 0.1 371/journal.pone.0095449.g002 



On day 7, over 70% of the bacterial community in sourdoughs 30- 
I, 30-11, and 30-III belonged to either L. hrevis, L. plantarum/ pentosus, 
or L. graminis/sakei/curvatus, respectively. However, on day 21 
fijrther changes in the composition of sourdoughs 30-1 and 30-III 
were detected. Over 90% of the bacterial amplicons identified in 
sourdough 30-1 belonged to L. paralimentarius /kimchii, while L. 
plantarum dominated in sourdoughs 30-11 and 30-111. On day 56, L. 
paralimentarius was the dominant species (>70%) in all sourdoughs 
fermented at 30°C, while the proportion of L. plantarum and L. 
brevis remained below 7%. In sourdough 30-1, L. crustorum was also 
found to comprise 15% of the bacterial community after 56 
backslopping cycles. 

Isolation, Identification, and Characterization of the 
Dominant Bacteria after 56 Backslopping Propagations 

In order to obtain pure cultures of bacteria after 56 days of 
propagation we randomly picked 1 0 colonies from MRS agar and 



10 from SDB agar for each of the six samples. In total, 120 
colonies were selected and analyzed by rep-PCR fingerprinting. 
Four different Rep-PCR fingerprint groups were detected (data 
not shown). Representatives of each group were identified using 
16S rRNA partial gene sequencing and were found to be L. 
plantarum [L. plantarum M30I-1, GenBank accession number 
KJ361844), L. brevis [L. brevis M30I-2, GenBank accession number 
KJ36I843), L. paralimentarius [L. paralimentarius M30I-3, Genbank 
accession number KJ361845) and L. crustorum [L. crustorum M30I-9, 
GenBank accession number KJ361846) with at least 98% identity. 
In most cases the fraction of each species identified in the 
sourdough samples are comparable with pyrosequencing data 
(Figure 4). 

Metabolic profiles of the four dominant LAB species were 
determined (Table 3). The strain L. plantarum M30I-1 was able to 
ferment the largest range of carbohydrates. Together with L. 
paralimentarius M30I-3, L. plantarum M30I-1 was able to ferment all 
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Table 2. Identification of DGGE bands obtained from 


six rye sourdoughs propagated at 20°C and 30°C during 56 days. 




Band number on 


Closest match in 






DGGE 


GenBank 


Similarity 


Accession Nr. 


1 


Gamma Proteobacteria species 


163/168 (97%) 


gb|GU352675.1 


2 


Pantoea agglomerans 


166/166 (100%) 


gb|KC355300.1 


3 


Pantoea agglomerans 


165/166 (99%) 


gb|KC355300.1 


4 


Pantoea agglomerans 


165/166 (99%) 


gb|KC355300.1 


5 


Enterococcus casseliflavus 


120/140 (86%) 


ref|NR_102793.1 


7 


Weissella confusa 


160/168 (95%) 


gb|KC845208.1 


8 


Lactobacillus graminis 


97/106 (92%) 


gb|KC836565.1 


9 


Lactobacillus curvatus 


133/155 (86%) 


gb|FJ609221.1 


10 


Lactobacillus curvatus 


143/155 (92%) 


gb|FJ609221.1 


11 


Pediococcus acidilactici 


1 1 8/1 24 (95%) 


gb|JF268323.1 


12 


Lactobacillus plantarum 


138/138 (100%) 


gb|JN863682.1 


13 


Lactobacillus curvatus 


247/250 (99%) 


gb|KF41 1435.1 


13 


Lactobacillus graminis 


247/250 (99%) 


gb|KF149819.1 


13 


Lactobacillus sal<ei 


247/250 (99%) 


gb|KF1 49680.1 


14 


Lactobacillus curvatus 


206/217 (95%) 


gb|KF41 1435.1 


14 


Lactobacillus graminis 


206/217 (95%) 


gb|KF149819.1 


14 


Lactobacillus sal<ei 


206/217 (95%) 


gb|KF1 49680.1 


15 


Pediococcus pentosaceus 


276/276 (100%) 


gb|JN851 781.1 


16 


Pediococcus pentosaceus 


254/257 (99%) 


gb|JN851 779.1 


17 


Lactobacillus sal<ei 


276/276 (100%) 


gb|JN851 763.1 


18 


Lactobacillus brevis 


257/257 (100%) 


gb|JN863690.1 


21 


Lactobacillus crustorum 


258/261 (99%) 


gb|KF1 93907.1 


23 


Lactobacillus brevis 


132/132 (100%) 


gb|KC845206.1 


24 


Lactobacillus plantarum 


140/140 (100%) 


gb|JN863682.1 


25 


Lactobacillus plantarum 


259/259 (100%) 


gb|KF31 8862.1 


26 


Secale cereale, complete genome 


121/122 (99%) 


gb|KC91 2691.1 


27 


Secale cereale, complete genome 


1 20/1 22 (99%) 


gb|KC91 2691.1 


28 


Lactobacillus paralimentarius 


113/113 (100%) 


gb|KC755102.1 


29 


Lactobacillus brevis 


129/129 (100%) 


gb|KC845206.1 


30 


Lactobacillus crustorum 


113/113 (100%) 


gb|KC755094.1 


31 


Lactobacillus paralimentarius 


114/114 (100%) 


gb|KC755 102.1 


doi:1 0.1 371 /journal.pone.0095449.t002 



four major cereal carbohydrates (glucose, fructose, maltose, 
saccharose) while L. brevis M30I-2 and L. crustorum M30I-9 were 
able to fermem either glucose and fructose or glucose, fructose, 
and maltose, respectively. 

Isolation, Identification, and Characterization of the 
Dominant Yeasts after 56 Backslopping Propagations 

After 56 days of propagation 10 colonies were picked from YPD 
agar plates for each sourdough sample (60 colonies in total) and 
fmgerprinted using M13 primer. Four chflFerent patterns were 
observed (data not shown) and representatives from each group 
were identified. They belong to Kazachstania mispora [K. mispora 
Y30I-10, GenBank accession number KJ361847), Candida glabrata 
(C. glabrata Y30II-1, GenBank accession number KJ361850), 
Saccharomyces cerevisiae (S. cerevisiae Y30II-9, GenBank accession 
number KJ361848) or Candida krusei {C. krusei Y30II-5, GenBank 
accession number KJ361849) based on 26S rRNA partial gene 
sequencing. K mispora was found to be the only dominant species 
in all sourdoughs fermented at 20°C, but was also found in 



sourdough 30-1 which was propagated at 30°C. C. glabrata was the 
only yeast species found in sourdough 30-III while sourdough 30- 
II contained three yeast species, S. cerevisiae, C. krusei, and C. 
glabrata, in the proportion 2:3:5. 

Isolated representatives of K. mispora, C. krusei, and C. glabrata 
were not able to ferment maltose and had a narrow carbohydrate 
fermentation profile compared with S. cerevisiae (Table 4). How- 
ever, only C. krusei and C. glabrata could ferment .A'-acetylglucosa- 
mine and trehalose, respectively. 

Discussion 

Previous research has established that stabilization of LAB 
consortia in spontaneously started sourdoughs occurs in a three- 
stage evolution process over the course of five to ten days [2,8,10- 
11,28]. During this time the acidity drops and stabilizes together 
with the LAB count to a level common in mature sourdough. The 
majority of these studies used liquid sourdoughs, which were 
fermented using continuous stirring. To our knowledge there are 
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Figure 3. Pyrosequencing of 16S rRNA gene amplicons from spontaneous rye sourdoughs propagated for 56 days. Three sourdoughs 
were fermented at 20 C (20-1, 20-11, 20-111) and three sourdoughs fermented at 30 C (30-1, 30-11, 30-111). Sourdoughs were sampled at days 1, 3, 5, 7, 21, 
56. The relative abundance at the species level based on partial 16S rRNA gene sequences is given. Species forming less than 5% of the population 
were grouped together and are shown as 'Others (<5%)'. (a) and (b) stand for two different sequencing runs. 
doi:1 0.1 371 /journal.pone.0095449.g003 
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Figure 4. Ratio of species in tlie sourdougKis fermented at 20 C (20-1, 20-11, 20-111) or 30'C (30-1, 30-11, 30-111) after 56 bacicsiopping 
cycles determined by plating on MRS and SDB media or by pyrosequencing of 16S rRNA gene amplicons. 
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Table 3. Carbohydrate fermentation 


profiles of Lactobacillus species isolated 


on day 56 of sourdougli bacl<lopping. 




Active ingredient 


Lactobacillus Lactobacillus 
plantarum brevis 
M30i-1 M30i-2 


Lactobacillus Lactobacillus 
paralimentarius crustorum 
M30i-3 IVI30i-9 


L-arabinose 


+ + 


+ - 


D-ribose 


+ + 


+ - 


D-xylose 


+ + 


+ - 


D-galactose + + — + 


D-glucose 


+ + 


+ + 


D-fructose 


+ + 


+ + 


D-mannose 


+ - 


+ + 


D-mannitol 


+ - 


- 


D-sorbitol 


+ - 


- 


N-acetylglucosamine 


+ + 


+ + 


Amygdalin 


+ - 


+ + 


Arbutin 


+ - 


+ - 


Esculin ferric citrate 


+ - 


+ - 


Salicin 


+ - 


+ + 


D-celiobinose 


+ - 


+ + 


D-maltose 


+ + 


+ + 


D-lactose (bovine origin) 


+ - 


- + 


D-melibiose 


+ + 




D-saccharose (sucrose) 


+ - 


+ - 


D-trehalose 


+ - 


+ + 


Inulin 


+ - 




D-melezitose 


+ - 


+ - 


D-raffinose 


+ - 




Gentiobinose 


+ - 


+ + 


Sodium Gluconate 


+ + 


+ - 


Sodium 5-Ketogluconate 


- + 





doi:l 0.1 371 /journal.pone.0095449.t003 



no studies that have monitored the fate of a bacterial community We followed 56 daily backslopping cycles of spontaneously 

after the sourdough has reached maturity under controlled started semi-solid rye sourdough at two temperatures with six 

laboratory conditions, where the raw flour is the sole bacterial parallels. While the sourdoughs achieved maturity in 10 to 12 
source. 



Table 4. Carbohydrate fermentation 


profiles of yeast species 


isolated on day 56 of 


sourdough bacl<lopping. 






Saccbaromyces 


Kazachstania 


Candida 


Candida 




cerevisiae 


unispora 


krusei 


glabrata 


Active ingredient 


Y30ii-9 


Y30I-10 


Y30il-5 


Y30il-1 


D-glucose + + + + 


Glycerol 






-1- 




D-galactose 


-1- 


-1- 






N-acetylglucosamine 






-1- 




D-maltose 


-1- 








D-saccharose (sucrose) 


-1- 








D-trehalose 








+ 


D-melezitose 


-1- 








D-raffinose 


-1- 








doi:l 0.1 371 /journa!.pone.0095449.t004 
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propagation cycles, further succession of LAB species was observed 
even after 42 cycles. 

The low number of sourdough specific LAB in rye flour could 
be the main reason for the observed instability. The concentration 
of bacteria in rve flour varies between 10* to 10^ CFU/e 
depending on the climate, time of harvest, and both miUing and 
storage conditions [29]. Using 16S rRNA gene pyrosequencing we 
found that the microbiota in the raw flour used in this study 
contained LAB species at a subdominant level [<!%) and the 
majority of these were species that typically do not dominate in 
sourdough communities. Thus the strains that dominated in the 
six mature sourdough paraUels may have been present in very low 
concentrations in the raw flour. The effect of introducing low 
numbers of sourdough competent LAB together with the small 
volumes of flour used to prepare the sourdoughs could together 
work to create an uneven distribution of bacteria between the 
three parallels conducted at the same temperature. As with our 
observations, Minervini et al. [31] observed the succession of LAB 
strains during propagation by comparatively following the 
microbial community of mature wheat sourdough during propa- 
gation in an artisan bakery and a controlled laboratory environ- 
ment. They attribute the cause of the observed succession to be 
differences between the batches of flour used. 

The temperature used for sourdough fermentation is one of the 
key factors determining the composition of the microbial 
community in sourdough [15,32]. The process of establishing 
sourdough microbial consortia occurred more rapidly at 30°C, as 
evidenced by the prevalence of LAB in these sourdoughs already 
after the first fermentation cycle. Several studies have shown that 
only one fermentation at 30°C is needed to enrich LAB in 
sourdough [11,32], whereas up to three renewal cycles were 
needed to detect LAB at a lower temperature (23°C) [32]. We also 
observed LAB after one fermentation at 2()°C, albeit at a 
significandy lower fraction of the total community compared with 
fermentations at 30°C. In addition, the normalized number of 
expected OTUs indicates that the species richness during the first 
three propagation cycles was higher in sourdoughs fermented at 
higher temperature. Interestingly, enterobacteria persisted over 
more fermentation cycles in bacterial communities of sourdough 
fermented at 30°C (up to three fermentations) than in sourdoughs 
fermented at 20°C. 

Facultative heterofermentative bacteria dominated in sour- 
doughs fermented at 30°C after 56 propagation cycles, while in 
sourdoughs fermented at 20°C both obligate and facultative 
heterofermentative LAB were dominant. It has been shown that 
despite the metabohc inefficiency of obligatory heterofermentative 
lactic acid bacteria, these species commonly dominate in 
sourdough fermentations [2]. Balance between homofermentative 
and heterofermentative lactic acid bacteria depends on dough 
yield, redox potential and fermentation temperature [30]. In this 
study fermentation temj)C'rature shifted the ratios between the 
limited number of LAB species in the community rather than 
select for different species. Rep-PCR analysis of isolates obtained 
from the sourdoughs after 56 renewal cycles revealed that 
fingerprints among isolates of the same species were similar. 

Pyrosequencing results establish that during the development of 
bacterial consortia in sourdough, irrespective of fermentation 
temperature, L. plantamm dominated or codominated in the 
majority of sourdoughs. However, after ,56 propagation cycles L. 
plantamm had been replaced by L. paralimentarius in most batches. L. 
plantamm is considered a highly acid-tolerant LAB that dominates 
in fermentation processes with vegetables and cereals due to its 
metabolic flexibility and low pH adaptation [33-35]. It also 
dominated in bacterial communities of four liquid spontaneously 



started laboratory rye sourdoughs after 10 backslopping renewals 
[8]. The carbohydrate fermentation pattern of L. paralimentarius 
isolated in this study is more constrained compared with L. 

plantamm, however, they are able to ferment all major carbohy- 
drates (maltose, saccharose, fructose, and glucose) [36]. On the 
other hand, succession within bacterial communities may depend 
on many other factors aside from carbohydrate metabolism, 
including amino acid metabolism and tolerance to acid stress 
conditions [10]. Several studies have shown that the competiveness 
of LAB often depends on their intraspecies diversity and is strain- 
specific [28,36]. 

Fermentation temperature also influenced the composition and 
diversity of the yeast community within the sourdoughs studied. 
Both yeast species detected in 30°C fermentations after 56 
propagation cycles, i.e. Saccharomyces cerevisiae and Candida krusd 
belong to the six most frequently encountered species in 
sourdoughs [37], while Candida glabrata is considered to be a 
prevailing species during liquid laboratory sourdough and teff 
fermentations [28,38]. Kazachstania unispora was the only yeast 
species found in sourdoughs fermented at 20°C. This species has 
been documented twice to exist in a sourdough ecosystem, albeit 
in low abundance: once in Belgian artisan wheat sourdough [37] 
and once in Finnish rye sourdough starter [39]. In contrast, K. 
unispora plays a significant role in both the ripening of cheese and 
in the production of fermented milk products such as a kefir and 
koumiss (reviewed by [40]). Additionally orange, sugarcane or 
mixed vegetable juices favour growth of K unispora. Ambient or 
low temperatures and high organic acid concentrations are specific 
for most of these processes [38-39,41-44]. No clear relationship 
was found between the yeast and LAB strains detected in this 
study. However, most isolated yeast species were maltose negative 
and trophic interactions between LAB and yeasts could be 
suspected because cooperation in utilization of maltose is the 
most frequently cited reason for the co-occurrence of yeasts and 
hAB (reviewed in [45]). 

Studies of microbial consortia in food matrices, such as 
sourdoughs, are commonly based on a combination of culture 
dependent and culture independent methods [46-47]. This study 
applies both approaches to analyze the bacterial community 
within mature sourdoughs after 56 propagation cycles. Plating on 
both SDB and MRS media, which are both commonly used for 
sourdough specific LAB isolation [48], provided similar numbers 
of colony forming units. In most samples species occurrence and 
ratio among identified isolates was comparable with pyrosequen- 
cing data. 

The results of both DGGE and 16S rRNA gene pyrosequencing 
were generally in agreement despite the fact that they amplified 
different regions of the 16S rRNA. The bacterial cUversity 
suggested by pyrosequencing results was greater than that revealed 
by DGGE because of the diflFerences in their detection limits. Both 
methods are based on DNA amplifK:ation and share similar 
limitations, e.g. possible inefficient DNA isolation from some 
organisms. However, high-throughput sequencing analysis is 
considered quantitative if efforts are made to minimize changes 
in the original proportion of microbial cells caused by DNA 
extraction [1 1]. Amplification of DNA from dead organisms could 
be critical for the first step of sourdough propagation which would 
overestimate the diversity of metabolically active bacteria in flour 
and sourdoughs. Similar distribution of OTUs between DNA and 
RNA samples has been shown for mature sourdoughs [11]. 

Both DGGE and pyrosequencing often fail to discriminate 
between closely related species due to both the insufficient length 
and accuracy of sequences. Current developments in high- 
throughput sequencing technologies are expected to overcome 
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these problems. Despite its limitations, high-throughput sequenc- 
ing has the potential to become a powerful tool for the culture- 
independent study of sourdough microflora since it offers a more 
in-depth analysis [49]. 

Supporting Information 

Table SI Characterization of pyrosequencing data obtained 

from the analysis of sourdough samples. 

(PDF) 

Table S2 Relative abundance (%) of partial 16S rRNA gene 
sequences obtained by pyrosequencing of spontaneous rye 
sourdoughs propagated for 56 days. 
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